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We study the information geometry and the entropic dynamics of a 3-D Gaussian statistical 
model. We then compare our analysis to that of a 2D Gaussian statistical model obtained from the 
higher-dimensional model via introduction of an additional information constraint that resembles 
" " ^ the quantum mechanical canonical minimum uncertainty relation. We show that the chaoticity 

I . (temporal complexity) of the 2D Gaussian statistical model, quantified by means of the Information 

(«~^ ' Geometric Entropy (IGE) and the Jacobi vector field intensity, is softened with respect to the 

I chaoticity of the 3D Gaussian statistical model. 
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; I. INTRODUCTION 

f~i I A very important problem in modern science concerns the description and the understanding of the elusive concept 
of complexity [ItlM- There are pragmatic reasons of primary importance in quantum information science that justify 
the study of complexity, for example the problem of quantifying how complex is quantum motion Q. Unfortunately, 
our knowledge of the basic connections between complexity, dynamical stability, and chaoticity in a fully quantum 
domain is not satisfactory 0, The concept of complexity is very difficult to define, its origin is not fully understood 
P-[Tl| and it is mainly for this reason that several quantitative measures of complexity have appeared in the scientific 
literature [H-O- In classical physics, complexity measures are settled in a much better way. The Kolmogorov- Sinai 
metric entropy (l2| , that is the sum of all positive Lyapunov exponents [isj , is a powerful indicator of unpredictability in 
^ ■ classical systems and it measures the algorithmic complexity of classical trajectories [l^ - [l7| . Other known measures 
Tljj- [ of complexity are the logical depth 0], the thermodynamic depth [l9|, the computational complexity [l^l, the 
stochastic complexity [2]| and many more. Ideally, a good definition of complexity should be mathematically rigorous 
t ■ and intuitive at the same time so that we are able to tackle complexity-related problems in computation theory and 
statistical physics as well. Of course, a quantitative measure of complexity is truly useful if its range of applicability 
is not limited to few unrealistic applications. For this reason, in order to properly define complexity measures, the 
purpose for defining such a measure and what it is intended to capture should be clearly stated. 
] ' Another problem of great theoretical interest is understanding how to compare quantum and classical chaos (tem- 
poral complexity) and explaining the reason why the former is weaker than the latter [23 - |23 |. Indeed, it is comm only 
conjectured that the weakness of quantum chaos may be a consequence of the Heisenberg uncertainty relation [22|-[24{. 
It is also known that a quantum description of chaos is qualitatively different from a classical description and that the 
later cannot simply be considered an approximation of the former. Indeed, the only trace of quantum theory which a 
■ classical description may retain is the canonical Heisenberg's uncertainty relation, namely a minimum spread of order 

fi" in the 2n-dimensional phase space [1^ {h '= ^ where h is Planck's constant). 

Motivated by such considerations and following the line of reasoning presented in ^6|, we will discuss in this 
paper the manner in which the degree of complexity changes for a statistical model (the probabilistic description 
of a physical system) in the presence of incomplete knowledge ("softening effects") when the entropic dynamics (or, 
information-constrained dynamics) [27j on the underlying curved statistical manifolds becomes more constrained. By 
" softening effects" we mean any attenuation in the asymptotic temporal growth of complexity indicators of motion. 
Our aim is to reduce the probabilistic description of dynamical systems in the presence of partial knowledge to 
information geometry (Riemannian geometry applied to probability theory, see '28']) and inductive inference [29 "si?]. 
To achieve such a task, we have developed a theoretical framework termed IGAC, Information Geometric Approach 
to Chaos (s^ . [35} , where information geometric techniques are combined with Maximum relative Entropy methods 
[29l - l33j to study the complexity of informational geodesic flows on curved statistical manifolds (statistical models) 
underlying the probabilistic description of physical systems in the presence of incomplete information. IGAC is the 
information geometric analogue of conventional geometrodynamical approaches to chaos [sB, [s^] where the classical 
configuration space is being replaced by a statistical manifold with the additional possibility of considering chaotic 
dynamics arising from non conformally flat metrics (the Jacobi metric is always conformally flat, instead). It is 
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an information-geometric extension of the Jacobi geometrodynamics (the geometrization of a Hamiltonian system 
by transforming it to a geodesic flow [38|). The reformulation of dynamics in terms of a geodesic problem allows 
the application of a wide range of well-known geometrical techniques in the investigation of the solution space and 
properties of the equation of motion. The power of the Jacobi reformulation is that all of the dynamical information 
is collected into a single geometric object in which all the available manifest symmetries are retained- the manifold on 
which geodesic flow is induced. For example, integrability of the system is connected with existence of Killing vectors 
and tensors on this manifold. The sensitive dependence of trajectories on initial conditions, which is a key ingredient 
of chaos, can be investigated from the equation of geodesic deviation (the so-called Jacobi-Levi-Civita equation or 
JLC equation) . Within the IGAC formalism, conventional signatures of the chaoticity of classical dynamics emerge 
via the asymptotic exponential divergence of the JLC vector field intensity and/or via the specific behavior of the 
asymptotic growth of the Information Geometric Entropy (IGE) (i.e., the logarithm geodesic volume) [26j . 

In view of the above-mentioned considerations, we discuss here the information geometry and the information- 
constrained dynamics of a 3D Gaussian statistical model. We then compare our analysis to that of a 2D Gaussian 
statistical model obtained from the higher-dimensional model via introduction of an additional information constraint 
that resembles the quantum mechanical canonical minimum uncertainty relation. We show that the chaoticity (tem- 
poral complexity) of the 2D Gaussian statistical model, quantified by means of the IGE [26] and the Jacobi vector 
field intensity, is softened (that is, attenuated) with respect to the chaoticity of the 3-D Gaussian statistical model. 
In view of the similarity between the information constraint on the variances and the phase-space coarse-graining 
imposed by the Heisenberg uncertainty relations, we suggest that our work provides a possible way of explaining the 
phenomenon of suppression of classical chaos operated by quantization. 

The layout of this article is as follows. In Section II, we present the basic differential geometric properties of both 
the 3D and 2D Gaussian statistical models. In Section III, we describe the geodesic paths on the curved statistical 
manifolds underlying the entropic dynamics of the three and two-dimensional statistical models. In Section IV, we 
study the chaotic properties of the information-constrained dynamics on the underlying curved statistical manifolds 
by means of the IGE and the Jacobi vector field intensity. Our final remarks appear in Section V. 



II. GEOMETRY OF THE STATISTICAL MODELS 



The statistical models studied are a 3D uncorrelated Gaussian model and a 2D Gaussian statistical model obtained 
from the higher-dimensional model via introduction of an additional information constraint that resembles the canon- 
ical minimum uncertainty relation in quantum theory. For a brief and recent overview on the IGAC, we refer to [26| . 
Note that the dimensionality {2D, 3D) pertains to the macroscopic variables. 



A. The 3D Statistical Model 



The probability distributions p (x, y\fJ-x, <^x, Cj,) that characterize the 31? Gaussian statistical model are given by. 



p{x, y\nx, Ox, cfy) 



dcf 



27r crj;Cr.y 



exp 



1 . ^2 1 2 



(1) 



with (Jx and Oy in and /ij, in M. The Gaussian here is two dimensional in its microscopic space (x, y) but 
three dimensional in its macroscopic (contextual or conditionally given parameters) space. We will now relax this 
conditionality to explore the space of Gaussians described by [i^ , c^. and Oy . The infinitesimal Fisher- Rao line element 



dsg^ reads, 

3 

I, m=l ^ ^ V 

where the Fisher- Rao information metric gim {0) is defined as psj . 



d\ogp{x, y\fix, Oy) d\ogp{x, o^, Oy) 



with 0= {e\ d\ 03) = inx, a-x, Oy) and where (6) has the following matrix representation. 



(3D) 



Z:^ \ 



9lm i^) - 















(2) 



(3) 



(4) 
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Using it follows that the non- vanishing Christoffel connection coefficients [39j . 



= ^5''™ + 9j5.m - d„,gij) , (5) 

are given by, 

pi - pi L - — - — - — (fw 

^12—^21— J i 11 — r, I 22 — 1^33— ■ IDJ 

The scalar curvature VP'^^ of the statistical manifold of probability distributions in ^ is given by, 

7^(3C) = (0) i?n + (^) + ^33 ^^^^ (7) 

with g^"'gmk = S^. and the Ricci curvature tensor Rij defined as [39| . 

„ dcf „ T^fe Q j^k I T^/c T^n T^mr^/c ^o^ 

Substituting © into (|S|), it turns out that the non vanishing components of Rij are, 

-^11 = -^22 = J- (9) 



Riemann curvature tensor i?"^p [3; 



Thus, the scalar curvature becomes VS^^^ = —1. We finally point out the only non- vanishing component of the 

-R/Ji/p = '^ij^'^p ^ '^pTJ^i/ + r^j^r^p — rj^pF^j^, (lo) 
R\i2 = -\- (11) 



is given by, 



As a final remark, note that 



in agreement with Eq. ([7|). 



7^(3^) = (i?i2i2 + i?2i2i) W g'' (0) = -1, (12) 



B. The 2D Statistical Model 

The probability distributions p (x, y; /i^,, cr) that characterize the 2D Gaussian statistical model are given by, 

(13) 



, . dcf 1 

p{x, y; fix, 0-) = TT^exp 



1 , ^2 2 



with (T in Rq and /ij, in K. The probability distribution p (a;, y; fix, c) may be obtained from p {x, ylux, <^x, o'y) with 
the addition of the following macroscopic constraint, 

UxOy = S^, (14) 

where Y? is a constant belonging to and Ox = cr. The macroscopic constraint (jl4p resembles the quantum 
mechanical canonical minimum uncertainty relation where x denotes the position of a particle and y its conjugate 
momentum. The infinitesimal Fisher-Rao line element rfs|^ reads, 

2 ^ . 

dslo^ E 5[™''^(^)rf^'rf^™ = ^rfM^ + ^rf^', (15) 

m— 1 
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where the Fisher-Rao information metric gim {0) is defined as, 

with = (0^, 6^) =^ (/ia;, cr) and where g'lf-^f ^ (0) has the following matrix representation, 

9lT m-^HD- (17) 



Using (|17p . it follows that the non- vanishing ChristofFel connection coefficients T^^ are given by, 

^12 = ^21 = 1 r?! = I r22 = — ■ (18) 

The scalar curvature HS^^^ of the probability distributions in ([T3l) is given by, 

7^(2^) = .gll (0) R,, + .g22 (0) i?^^ = (19) 



with g Qrak = (^fc and where the only non- vanishing Ricci curvature tensor component Rij is, 

1 fi. _ 1 

Observe that the only non- vanishing component of the Riemann curvature tensor is 



^11 — —-rZ2' ^22 — 2' (^^) 



^212 = 2 ' (21) 



7^(2^' = (i?i2i2 + i?2i2i) .9" (0) g'' {9) = -i, (22) 



thus. 



in agreement with Eq. (|19p . From ([T2|) and ([22|) . it turns out that the 3D statistical model is globally more negatively 
curved than the 2D statistical model. 

III. GEODESIC MOTION ON CURVED STATISTICAL MANIFOLDS 

In this Section, we present the geodesic paths on the curved statistical manifolds underlying the entropic dynamics 
of the three and two-dimensional Gaussian statistical models. Such paths are obtained by integrating the geodesic 
equations given by {33] , 



d 9'' (s) d6^ d0 



dT^ W 

Ira 



+ rL(^)-7--7- = o, (23) 



where F^^ {9) are the Christoffel symbols defined in ([5]) 



A. The Three-dimensional Case 

Substituting ^ into ([25]) . the set of nonlinear and coupled ordinary differential equations in (|23p reads, 

Q _ <PjJ^ 2 d^x dax 
dr^ Ox dr dr ' 

^ _ d'^Gx ^ J_ f dfJ.x \^ 1 f da 



dr^ 2ax \ dr J Ux \dT 



0=^--f^V- (24) 
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A suitable family of geodesic paths fulfilling the geodesic equations above is given by (technical details appear in 
Appendix A), 



(/io + 2a-o) [l + exp (2croA+T)l - 4cto 2cro exp (ctoA'^^t) 



1 + exp (2(ToA'^t) 



1 + exp (2croA+r) 



and, 



c^yM = o-oCxp(-A/T), 



(25) 



(26) 



where /io '= fJ,x (0), ctq =' f^x (0), CTq =' ay (0), e (r/), A/ =' log (^) and A'+ G M+ (see Appendix A) 



dcf 



dof 



dcf 



dcf 1 



B. The Two-Dimensional Case 



Substituting (|T8|) into ([23l) . the set of nonlinear and coupled ordinary differential equations in (|23|) reads, 







dr^ a dr dr' 







1 /daV 
a [d^ 



A suitable family of geodesic paths fulfilling the geodesic equations above is given by, 

(^0 + 2a-o) [1 + exp (2a-oA+T)] - Acfq 



fJ-x (t) = 



and. 



a(r) = 



1 + exp (2(ToA+t) 

2^0 exp (ctqA+t) 
1 + exp (2(ToA+t) ' 



(27) 



(28) 



(29) 



where fig =^ fix (0), ctq =^ cr (0) and A+ = ^ belongs to M+ (see Appendix A). 



IV. INDICATORS OF CHAOTICITY 

In this Section, the chaotic properties of the information-constrained (entropic) dynamics on the underlying curved 
statistical manifolds are quantified by means of the Information Geometric Entropy and the Jacobi vector field 
intensity. The relevance of such quantities as suitable indicators of chaoticity was also investigated in [i^l ■ 



A. Information Geometric Entropy 



We point out that a suitable indicator of temporal complexity within the IGAC framework is provided by the IGE 
Sms (r) Si, 



Sm, (t) =^ log vol 



The average dynamical statistical volume vol 



V. 



(geodesic) 

e 



(t) 



vol 



V, 



(geodesic) 



ir) 



'^'^^ lim ( — I dr' vol 



r— >oo \ T 



p(gcodcsic) 



is defined as, 

(geodesic) ^^/-j 



Pi 



(30) 



(31) 
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where the "tilde" symbol denotes the operation of temporal average. For the sake of clarity, we point out that in the 
RHS of (|3ip . we intend to preserve the temporal-dependence by considering the asymptotic leading term in the limit 



of r approaching infinity. For a n-dimensional manifold, the volume vol 



(geodesic) , i\ 



(r') in ([3T|) is given by. 



vol 



^-^(geodesic) ^^i^ 



def 



=-'(r') 



(32) 



where P(Ms, g) (^^ v: is the so-called Fisher density and equals the square root of the determinant of the metric 
tensor g^j,^ {9), 



P(Ms, g) 



) ^g{[e\...,e-)). 



(33) 



The integration space T)^ 



(geodesic) , r\ ■ 



(r') in (|32l) is defined as follows. 



(geodesic) / /\ dcf 



(r') = {e = {e\..., 0") : (0) <e'' <e'' (r')} , 



where k = 1,.., n and 6^ = 9'^ (s) with < s < r' such that, 



0. 



(34) 



(35) 



ds^ ds ds 

(geodesic) (jg^ is a n-dimensional subspace of the whole (permitted) parameter space 



The integration space 

Vf°'\ The elements of 15^8°°''^^'^' (r') are the n-dimcnsional macrovariables {9} whose components 9'^ are bounded by 
specified limits of integration 9'^ (0) and 9'^ (r') with k = 1,.., n. The limits of integration are obtained via integration 
of the n-dimensional set of coupled nonlinear second order ordinary differential equations characterizing the geodesic 
equations. Formally, the IGE Sm, {t) is defined in terms of a averaged parametric {n + l)-fold integral (r is the 
parameter) over the multidimensional geodesic paths connecting 9(0) to 9{t). 



In the cases being investigated, using and it turns out that the Sj^ 
reads, 



313 Model) 



of the SD statistical model 



S 



{3D Model) 
Ms 



1 , ,(3D Model) 



with. 



V 



(3D Model) 
Ms 



(r) 



1 exp (ctqA^t) 



(2(70 + /io) CToA/A'+T + (2ao - Mo) aoXfX'+Te-^""^'^- 
- (A/ + A^o-Q In a'g) (2(To + Mo) + 
- ((To a; Ina^ - A/) (2ao - Mo) e-^-orA' 



In the asymptotic limit, we get 



V 



(3D Model) f n -rgjl 



Ms 



that is. 



A/ Mo + 2(70 \ / \/ \ 
-f • -2 exp (aoA+r) 



(36) 



(37) 



(38) 



c,(3D Model) / s t»1 , , 



Similarly, using ([281) and ((M]), it follows that the S 



(2D Model) 



Ma 



of the 2D model becomes. 



S 



(2D Model) 
Ms 



, , ,(2D Model) / n 



with. 



(2D Model) 



Ms 



ir) 



1 (mo + 2ao) -I- (2ao - ^p) exp (-2goA+r) 
A+ctq r exp (-(TqA+t) 



(39) 



(40) 



(41) 
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In the asymptotic limit, we obtain 



V 



(2D Model) 
Ms 



ir) 



T>1 



/io + 2cro~\ exp((7oA+r) 



(42) 



1 o(2-D Model) , 

and bj^ ' becomes, 



(,(2D Model) 



ir) 



r»l 



A+T. 



Combining (|43| and (|39|) . it finally turns out that 



o(2D Model) / n t»1 



•5 



(3D Model) 
Ms 



ir) 



with -— t = < 1. 

a; V2 



(43) 



(44) 



Eq. pi)) is quite interesting since it quantitatively shows that the IGE is softened when approaching the two- 
dimensional case from the three-dimensional case via the introduction of the macroscopic constraint ((Ti)) that is 
reminiscent of Heisenberg's minimum uncertainty relation where x denotes the position of a particle and y its conjugate 
momentum. 



B. Jacobi Vector Field Intensity 

The Jacobi-Levi-Civita ( JLC) equation of geodesic deviation is a complicated second-order system of linear ordinary 
differential equations. It describes the geodesic spread on curved manifolds of a pair of nearby freely falling particles 

travelling on trajectories 9'' (r) and 6'^ (r) =^ 6'' (r) + 66'' (t). The JLC equation is given by (42j . 



with k = 1, 2, 3 (in the three-dimensional case) and where the covariant derivative ^^j^l^^ along the curve 6"^ (r) is 
defined as, 

The Jacobi vector field components J*^ are given by, 

J^^^A.^^^=^' (^^^^) SX„ (47) 

and Rap-yS is the Riemann curvature tensor defined as 42] , 

Kup d.Kp - ^P^M. + - ^pKu ■ (48) 

The covariant derivative in (|45|) is defined as (43l |. 

-Dr^-^^ ^^"^^^ + ^"^^^17^ + ^"^^ ^^^-^ W.^^^^- (49) 

Equation (j45p forms a system of coupled ordinary differential equations linear in the components of the deviation 
vector field (j47p but nonlinear in derivatives of the metric tensor {9). It describes the linearized geodesic flow: 
the linearization ignores the relative velocity of the geodesies. When the geodesies are neighboring but their relative 
velocity is arbitrary, the corresponding geodesic deviation equation is the so-called generalized Jacobi equation [il j. 
The nonlinearity is due to the existence of velocity-dependent terms in the system. Neighboring geodesies accelerate 
relative to each other with a rate directly measured by the curvature tensor Rap-yS- 

Omitting technical details (see Appendix B) and setting A^d XsdO'o with Asd = A'^, the JLC equation for 
becomes, 

2A3i3-^ - 8A3D exp (-Aacr) — - 8A3D exp i-2A3nT) - SA^^ exp (-Agcr) = 0. (50) 



ar"^ dr dr 
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Similarly, it can be shown that the JLC equations for and read, 



dp 



+ 4A3D exp (-Aacr) — + 2A3U— + KioJ = 0, 



(51) 



and, 



+ 2A/- 



(52) 



dr2 ' dr 

respectively. In the asymptotic limit, it can be shown that the three equations to integrate become (see Appendix B) 



dr2 



2A3D^+A2^J2 =0, 



2X 



dJ^ 
~d7 



0. 



(53) 



The asymptotic solutions are given by. 



(r)=cf^+C^'^exp(-2A3DT), 



J' (r) = ) exp {-Kwt) + C^'Vexp (-As^r) , 



(r) = Cf ^ exp (-A/r) + C^^'Vexp (-A/r) , 



(54) 



where Cy^ with fc = 1, 2, 3 and fc' = 1, 2 are real integration constants. The Jacobi vector field intensity J]^' (t) is 
defined as. 



t(3£i) 



J^f ) (r) 



[jHr)]\2[jHT)Y , 2[jHr 



alir) aUr) oUr) 



(55) 



Eqs. (Uni), (Uni and ^ imply that, 



(1) 



4a2 



■exp 2A3DT , ^ , ; 

f^l (t) 



(2) 



2a2 



-T^ and. 



C. 



(3) 



(56) 



It finally follows that J^^"* (t) reads. 



(r)«^exp(A3Dr) 



(57) 



Consider now the two-dimensional statistical model. Omitting technical details (see Appendix B) and setting K2D 
A2DC0 with \2D = A+, the two JLC equations for and read, 



dof 



rf2jl (ijl rfj2 

+ 2A2D— - 8A2i> exp (-A2DT) — - 4A2^ exp (-2A2dt) - SA^^ exp (-Aa^r) = 0, (58) 



and, 



-— + 2A2D exp (-A2DT) — + 2A2D— + A^^ J2 = 0, 
dr^ dr dr 



(59) 
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respectively. Following the line of reasoning provided for the three-dimensional case, the asymptotic integration of 
([55)1 and ^ lead to 



Ji (r) ^ C['^ + C^'^ exp (-2A2cr) 



(1), 



and, 



(r) = ^ exp (-Asdt) + C^'Vexp (-A2Z5T 



(2), 



(60) 



(61) 



respectively where C^., with k ~ 1, 2 and fc = 1, 2 are reaZ integration constants. The Jacobi vector field intensity 
Jj^^ (r) is defined as, 



Eqs. ([291), ([60]) and (|6T]) imply that. 



[jHr)]' 4[j^(r)] 
a2 (r) (72 (r) 



(62) 



^2 (r) 



(1) 



■ exp (2A2_Dr) 



4[J2(7 



a2(r) 



ft then follows that. 



J' 



{2D) 



Ms 



(1) 



2ao 



exp (A2DT) 



Thus, we have shown that in the asymptotic limit. 



J 



(3D) 



Ms 



(1) 



c. 



(2) 



■ exp (Asdt) and, J_^^^ (r) 



2crr 



that is. 



J 



(2D) 



-(A3D-A2D)r , t(3D) 



2crr 



■ exp (A2DT) , 



^vM? (^) with A3D - A2D > 0. 



(63) 



(64) 



(65) 



(66) 



Eq. ()66p is quite enlightening since it shows that the Jacobi vector field intensity is softened when approaching the 
two-dimensional case from the three-dimensional case via the introduction of the quantum-like macroscopic constraint 

m- 

We emphasize that our main findings (see Eqs. (j44p and (|66p ) presented in this work are limited to the asymptotic 
behavior at infinity on the r-scale of the selected indicators of cliaoticity, that is, the IGE and the JLC vector field 
intensity. However, recall that in quantum chaos, the shortest random time scale (or Ehrenfest time scale) tr is 
approximately given by (45|, 



(67) 



where A is the maximum Lyapunov exponent of the system, ^ =^ ^ is the Dirac constant, h is the Planck constant 
and 5* ~ J pdq is a characteristic reference value of a classical action. This time scale is especially important because 
on time scales of this order, the complete transition from quantum to classical dynamical chaos is allowed. Stated 
otherwise, on time scales of this order, quantum motion is similar to the classical one including the exponential 
instability . Therefore, it would be worthwhile investigating whether or not for the statistical dynamical Gaussian 
models studied here, there is any information geometric analogue of the standard random time scale, say Tr, such 
that one may find softening effects on a r-scale longer than a finite r^. This investigation would present two delicate 
points: first, we would need the JLC-equation analysis on finite-r scales which would be extremely difficult from 
a computational point of view; second, our geodesic affine parameter r is not the conventional time t. Therefore, 
unless the connection between r and t is clearly specified, this analogy would be admittedly vague. However, we 
believe that this investigation may be successfully tackled within the IGAC framework at least in specific cases, for 
instance for conservative chaotic Hamiltonian systems. In such cases, there is a neat connection between the standard 
time-scale and the geodesic parameter-scale and the difficulties in the integration of the JLC equation may not be 
insurmountable (47| . We will examine this issue in forthcoming works. 
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V. CONCLUDING REMARKS 



In this article, we studied both the information geometry and the entropic dynamics of a 31? Gaussian statistical 
model. We then compared our analysis to that of a 2D Gaussian statistical model obtained from the higher-dimensional 
model via introduction of an additional information constraint that resembles the quantum mechanical canonical 
minimum uncertainty relation. We showed that the chaoticity (temporal complexity) of the 2D Gaussian statistical 
model, quantified by means of the Information Geometric Entropy (IGE) and the Jacobi vector field intensity, is 
softened with respect to the chaoticity of the 31? Gaussian statistical model. Specifically, Eq. (|44|) shows that the 
IGE is softened when approaching the two-dimensional case from the three-dimensional case via the introduction 
of the macroscopic constraint (jl4p that resembles the quantum mechanical canonical minimum uncertainty relation. 
Furthermore, Eq. ()66|) confirms that also the Jacobi vector field intensity is softened when approaching the two- 
dimensional case from the three-dimensional case via the introduction of the macroscopic constraint ([T4l) . 

We stress that our information geometric analysis could be further generalized to accommodate non minimum 
uncertainty-like relations. However, such extension requires a more delicate analysis where Maximum relative Entropy 
methods are used to process information in the presence of inequality constraints (48j . Our work is especially relevant 
for the quantification of soft chaos effects in entropic dynamical models used to describe actual physical systems when 
only incomplete knowledge about them is available [49j. Furthermore, although we are aware that our analysis is 
not manifestly "quantum", our findings lead us to support the conjecture that quantum chaos is ultimately weaker 
than classical chaos because of Heisenberg's uncertainty relation, the most important difference between classical and 
quantum physics. Of course, a deeper analysis is needed and we leave it for future investigations. 
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Appendix A: Integration of Geodesic Equations 

We show here the details leading to the geodesic paths presented in Section III. 
Consider the following set of coupled nonlinear differential equations. 

Pa. a; _ Q ^^^^ _x ^ J-x_ -7^ =0. (Al) 



dr^ Ux dr dr ' dr^ 2(Jx \ dr J Ux \ dr 

Setting pLx =^ and dx =^ ''"'^dT^' i the two equations in ()Aip become. 



From (IA2I) we observe that 



■jlx - 2— fix = and, ax + tt^A' - — = 0- (A2) 

Cx ^Ox Ox 



^ = 2^ . (A3) 

Therefore, we have that 

fix (r) - A.al (r) , (A4) 
with Al E K. Substituting (jA4|) in the second equation in (|A2p . we get 

axox -<jI + = 0. (A5) 

Our goal is to integrate Eq. (jA5l) , find Ox (r) and finally compute fix (t) using (jA4l) . For the sake of simplicity, let us 
put ax (t) = y (r) and a =^ 4^ G Mg . Then, Eq. (jASp becomes 

yy-f + ay^ = 0. (A6) 
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As a first change of variables, let us set 

Substituting (HZj into we obtain 

As a second change of variables, let us set 

Therefore, it follows that 



y (t) = =x{t). (A7) 



XX - x^ + ai* = 0. (A8) 



dx (r) dci , . , 
x = ,^ ' = z(x). (A9) 



i = zz\ X = (z"z + z, z' = (AlO) 
^ ' ax 

Substituting (jAlOl) and (|^ in we get 

z" + az = 0. (All) 

Integration of (lAlip yields, 

z {x) = Ci sin (V^x + C2) , (A12) 
where the integration constant coefHcients Ci and C2 belong to M. Recalling that x — z {x) , we get 

^ da;'= / dr'. (A13) 



Cism{y/Ex' + C2) 
Upon integration, (jAlSp becomes 

/ax + C2 



■log 



tan 



r + Cs, (A14) 



with C3 S R. Solving for x = x (r), we finally obtain 

a; (r) = ^ {2 arctan (exp [V^Ci (r + C3)] ) - C2} . (A15) 
V a 

Finally, (r) reads 

2ciexp(ci^T + ciC37Sj) 

(^) = 7-^ 77, r- I o 

1 + exp {2ci,/c^T + 2clCz^/o2) 

where we have set ci — Ci, C2 ~ a ^ and C3 = C3. It is straightforward to verify that indeed (r) in (IA16[) 
satisfies the nonlinear differential equation (IA6|) . Finally, we can compute ji^ (t) using (jA4l) . It follows that, 

, \ nr r 2 , l^ I C4 [l + exp(2civ^T + 2ciC3Vc2)] -4ci 

/X, (r) = V4a / (r') dr' = ' , ^ ^ ^ „ ^ • (A17) 

J 1 + exp (2ci^/ci^r + 2ciC3^ci^j 

Therefore, the geodesic paths are given by 

C4 [1 + exp {2ci^T + 2ciC3V^)] ~ 4ci 2ci exp (ci ^^^^r + ciC3V^^) 

= T"^ 77 r- I o 1=\ " 7"^ 7= — T^^ 1=\ ■ ^^^^> 

1 + exp (2ci^/ci^r + 2ciC3^/ci^j 1 + exp (2ciy/c^T + 2ciC3^c^j 

As working hypothesis, we consider geodesic paths with C3 = and assume that the initial conditions given by 
fJ-x (0) = Ho and a (0) = ctq- The geodesies in (|A18|) become, 

(/io + 2cto) [1 + exp (2croA+r)l - 4cro , , 2cro exp (ctqA+t) 

/^^(■^) = 7-, 77 — TTT ^(^x{t) = — J- — — ^, (A19) 

1 + exp (2(ToA_^Tj 1 + exp {2(Jo'^+t) 
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where A'^ = ^Jcq, = \fa > since a =' ^ > 0. 

Using the same tricks (invertible changes of variables) presented above, it becomes straightforward integrating the 
last differential equation characterizing the three dimensional case, 



Its solution ay (r) is given by, 



<^y (t) = croexp(-A/T) , 



(A20) 
(A21) 



with ay (0) ^ a'o, ay (Tf) ^ e and A/ =^ 77 log ( ^ 



Finally, following the line of reasoning presented above and assuming the very same working hypothesis and initial 
conditions, it turns out that the geodesic paths in the two-dimensional case read. 



(Mo + 2f7o) [1 + exp (2(ToA+t)] - Aao 

Ma; (T) = : 7Z : ^ , a{T] 



1 + exp (2o'oA+r) 



2(70 exp (q-qA+t) 
1 + exp (2(ToA+t) ' 



where A+ = \/a > since a =^ ^ > 



(A22) 



Appendix B: Integration of Jacobi-Levi-Civita Equations 



We present here few technical details relative to the JLC equation analysis introduced in Section IV. 
The JLC equations in the three-dimensional case are three. The JLC equation for reads. 



Dr2 



ji I dax \ _ j2 d^J^x dax 



dr 



where we have used the relation i?2i2 = Recall that 



For we get. 



dJ^de^.^, jo^d^ 
dr dT ^ dT^ 



dr dr 



= 0, 



dO" de^ 



P) r^" 7^" . TP 



1 fdj^de^ dj^de^ 

dr^ ' ^^[~d^~d^ ^ ~d7~d^ 



r]Jj 



■2/32 



dT^ 



dT dT 



dT^ 



P'' dT dT 



(Bl) 



(B2) 



d0^ 



d9^ d9^ 



dT J dT dT 

where we have used the identity, 

dd"^ 



J' 



1 .2 fde 



dT 



+2(1 i2j +I12I11 1 



de'^ 



since Tjj = ^21 = r22- From (IBip and (|B3|) and after some algebra, we get 







2pi dax dJ^ ^ d^x dJ^ ^ 
dT^ dT dT dT dT 



J\ (B3) 



(B4) 



1: 



(r}2)^ 



dr 



dax 



+ 



J- 12 



o T^i d^J,x dax , „ /T^i \2 d^x dax 



dr dr 



(B5) 
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that is, 



dr2 dr dr dr dr 



J' 



Observing that, 



Eq. (IB6P becomes, 



12 j ^212 J 



(ry' - ^, ^212 = and, d^rl, + 2 (r}^)' - rI,, = 4' 



dr2 dr dr dr dr 



(B6) 



(B7) 



ri2^ + a2ri2(^ 



d/i^ 
dr 



that is. 



J' 



^1 cP/^ 4 d/i^dCT^ 
12 dr2 a2 dr dr 



dVi 

dr2 



2 dcrj;\ dji 
(7t dr / dr 



2 d^^\ dj2 
ax dr y dr 



(B8) 



dr2 



dcTj. 
"d7 



1 

2al 



dr 



J' 



1 d^^x 4 d^x dffx 



Ux dr^ Gx dr dr 

At this point we recall that the geodesic paths ^x [t] and Ox (r) for the iD statistical model are given by, 

(^0 + 2cro) [1 + exp {2^/a(7oT)] ~ 4ao 



fJ-x (r) 



and, 



o-x (t) 



1 + exp (2v^(7or) 

2ao exp (yacror) 
1 + exp {2^/aaoT) 



(B9) 



(BIO) 



(Bll) 
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From (|B16|) and (|B17|) . the JLC equation for becomes after some algebra, 
dr^ V dr dr dr dr 



dr dr ^"^ dr dr 



r?i + d2Tli— r i Hi 12 



pi 

^22^11"; 1~ -^^121 

dr dr 



dd^ dd^ 
dr dr 



d^0^ 



It 



o 1 / dO^ 



Recalling that 6^ — ji^ and 6*^ ^ Cx, we get 



dV2 
dT2 



2r 



11' 



dr dr 



+ 2T 



22' 



dax dJ^ 
dr dr 



2 ^^02 

dr 



dO^ 

d7 



2 d^ fix 



J2 



(Po-x 



Noting that the following relations hold 

92^11 + r^r^j + r22r2j 



(a2r2, + (r2,)^) + {rj^rl, + rI,,) 



R 



121 



-, i92r22 + (r22; 



2 \2 _ 2 



and, Tl^T\^ + i?22^ 



Eq. (|B19P becomes, 



c!2j2 / I dfix\ dJ^ f 2 d(Jx\ dJ 



dr"^ \<Jx dr I dr 



(J r dr J dr 



1 d^ fix 1 d^,x dox 



2a-r dr2 a'i. dr dr 



From (|B12|) we get. 



and, 



+ J2 



1 (i2cr;^ 2 f d(Tx\^ 1 / c^A'x ^ ^ 



(Tj, c?t2 0-2 \^ dr 



dr 



— ^7^ « 4:y/aao exp (-Vao-QT) , ^ 2y/aao, 

ax dr (Tx dr 



1 d'^fj.x (r) 1 d/Ua; (r) dcTj; (t) 



2(Ta. (ir2 a?. dr dr 



1 rf2o-^ (^t) 2 f dax {t)\ If dfix (t) 



ax dr^ ax \ dr 
Substituting (|B22)) and ((B23| into ((B2T|) . we obtain 

d2j2 



dT2 



d7 



AAsD exp (-A315T) ^ + 2A3i5^ + A2^ j2 ^ 



dj2 
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Finally, the JLC equation for is given by, 



Notice that, 



0. 



(B25) 



Therefore, the JLC equation for becomes 



^•^ dr y dr 



1^ 



(B26) 



dr2 



+ 



2 rfj3 ^ ^3 



(Ty rfr / dr 



1 (i^CTj, 2 / ddy 



Gy dr"^ Gy \ dr 



having recalled that. 



r^g = - — and, ^rig + (rig)' = \. 



(B27) 



(B28) 



Observing that, 



dT 



-a'X 



dr 



(To Aje ^ 



C^'^^i/ (^) 'i2„-rAf 2 d(Jy[T) 

^CToAfe ' 



the JLC equation for becomes 



ay dr 



1 dVjr) , 2 ( dGy[T) 
a„ + al \ dr 



(B29) 



(B30) 



Following the same line of reasoning, we are capable of computing (and integrating) also the JLC equations for the 
2D statistical model. 
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